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Abstract 

We study in this paper the properties of a gas of fermions interacting via a 
scalar potential v{q) = A-kc^ /q^ for g < A << kp at dimensions larger than 
one, where A is a high momentum cutoff and kp ''\s the fermi wave vector. In 
particular, we shall consider the ^ oo limit where the potential becomes 
confining. Within a bosonization approximation, effective Hamiltonians de- 
scribing the low energy physics of the system are constructed, where we show 
that the system can be described as a fermi liquid formed by chargeless quasi- 
particles which has vanishing wavefunction overlap with the bare fermions in 
the system. 
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I. INTRODUCTION 



In the last few years, there has been enormous interests in the study of U{1) gauge 
theories of fermionic systems in dimensions higher than one, as a result of appearance of 
effective gauge theories in the t — J type models and also in the studies of the = | 
Fractional Quantum Hall state 0,^. In particular, the behaviour of the systems in the 
confinement regime where charge excitations are " confined" are of special interests. In this 
paper we shall apply a bosonization procedure to study a system of charge e spinless fermions 
interacting with scalar gauge field in dimensions higher than one, with effective long 
distance gauge field action = (V0)^. The problem is equivalent to fermions interacting 
via scalar potential v{q) = Ane'^/q'^ for g < A, where A << kp is a high momentum cutoff, 
kp is the fermi wave vector. In particular we are interested at the ^ oo limit where 
the potential becomes confining. In this limit any non-uniform charge density fluctuations 
< p{q) >7^ with any wave vector g < A cost infinite electric field energy. As a result any 
physical state {ip > which survives in this limit satisfies the constraint piqjlip >= 0, where 
p{q) is the charge density operator. 

In this paper we shall show that within a bosonization approximation the projection to 
the physical states \ip > where density fluctuations are forbidden can be carried out in the 

— s> oo limit, the bosonization approximation being a generalization of usual bosonization 
method in one dimension to dimensions higher than one 0-0. Within the approximation 
effective low energy Hamiltonians which describes the particle-hole as well as single-particle 
excitations in the system can be constructed. In particular, the single particle excitations are 
described by operators which commutes with charge operator p{q), and has zero-overlap with 
the bare fermion operators in the system. The ground state of the system can be described as 
a fermi liquid of these new single particle operators and corresponding to a kind of 'marginal' 
fermi liquid in the original fermion description. The organization of our paper is as follows: 
in section two we shall outline our mathematical formulation of the bosonization procedure 
where the nature of the bosonization approximation will be explained. In section III we 
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shall study in detail the eigenstates and eigenvalues of the particle-hole excitation spectrum 
in the bosonization approximation for arbitrary values of coupling constant e^. In particular 
we shall show how the projection to the physical Hilbert space where charged excitations 
are forbidden is archieved in the — > oo limit. In section IV we shall consider the single 
particle excitations where 'physical' single particle operators which commute with density 
operators are constructed and the equations of motion for these single-particle operators 
are derived. We shall show that the system can be described as a marginal fermi liquid 
of the original fermions. Our results will be summarized in section V where some further 
comments will be given. 



II. MATHEMATICAL FORMULATION 

We consider a gas of spinless fermions interacting via scalar potential v{q) = Ane^ /(f in 
momentum space. The Hamiltonian of the system is, 

H-T.^nftk+lL E <q)p{q)p{-q), (1) 

where — {k)^^/2m and /(/''')^'s are fermion annihilation (creation) operators. We have 
set ^ = 1 to simplify notation. p{q) — Z^^/^^^2/fc-q/2 density operator for the 

fermions. L'^ is the volume of the system. A << fcj? is a high momentum cutoff. Notice 
that the Hartree {q — 0) interaction energy does not appear in the Hamiltonian as in usual 
Coulomb gas problem. The difference between our model and usual Coulomb gas problem 
has to be emphasized here. In usual Coulomb gas problem the high momentum cutoff A 
is taken to be infinity, or satisfies A » kp- In this limit a Wigner crystal is expected 
to be formed (in 3D) when becomes large because the potential energy term becomes 
dominating at length scale ~ kp^. In our problem where A << kp, the interaction is 
effective only at length scale >> inter-particle spacing and formation of Wigner crystal 
is not possible since the fermions are essentially non-interacting at distance scale of inter- 
particle spacing. Instead the system should remain in a liquid state and a new treatment 
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to the problem is necessary. We note that bosonization method is a natural tool to attack 
the problem in this limit The bosonization procedure can be formulate most easily 

by introducing the Wigner function operators p^(g) = f-^j^^i2fk-q/2- shall work in the 
path-integral formulation where the Wigner operators are introduced in the action of the 
system at imaginary time through Langrange multiplier fields, 



S= dr 
Jo 



+ 



d 



k+q/2 ^ 



g^O,q<A,k,k' 

where n is the chemical potential. A^(g) are Lagranger multiplier fields introduced to enforce 
the constraint that the Wigner operators are given by Pj^iq) = /^^/g'^fc-'f/s' particular, 
the original Hamiltonian (|l|) is recovered once the A^(g) field is integrated out. 

Integrating out first the fermion fields /(/^)'s we obtain an action in terms of Pj;{q) and 
A^(g) fields. 



S = —Fq — Trln 

P 



1-GoX 



^ E ^H^pM 

k,q^O,q<A 



1 



2L<^ 



E ^(?)Pfc(^)Pfc'(-^). (3) 



g^O,(J<A,fe,fc' 

where Fq is the free energy for an non-interacting Fermi gas. Gq and A are infinite matrices 
in wave vector and frequency space, with matrix elements given by 

1 



k,k 



, = Sk,k'9o{k) go{k) 



(4) 



and 



A 
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where k = {k, iuon) and = P- The Trln [1 — GqA] term can be expanded in a power 
series of i\j:{q) field. 
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GnA 



O(A^). 



Keeping terms to second order in A (Gaussian approximation), we obtain 



Trln 



l-GnA 



(6) 



Notice that the first order term in A gives the usual Hartree self-energy and is excluded in the 
problem. The «A^(g) fields in Action can be integrated out in Gaussian approximation, 
resulting in an quadratic action in terms of Pj:{q) fields only. We obtain 



E 



where 



k,k',q,iu>n 



1 



1 



{L%t.)+v{q) 



(7) 



r, E 9o{k + q/2, + i^n)9o{k - q/2, zfij = ''^^-^72 ''Um ^ 



iflr, 



= 0{—Ck) at zero temperature is the free fermion occupation number. Sp can be expressed 
in terms of canonical boson fields by introducing 



(9a) 



where Ar(g) = n 



k-q/2 ^k+q/2 



Correspondingly, we also have 



Ag(g-)| {e{A^,{q))af^{q,zu^) + e{-A^,{q})ai{-q, 



(9b) 



Putting eqs. (pf) back into Sp, we obtain after some straightforward manipulations. 



•^p = ^ E i-'^^n+^-^)a^{q,iuJn)af:{q,iuJn) + ^ 

k,q,iiL>n k,k',q,iLU„ 



y: v{q)MM^k'm (10) 



+ a-fc(-9, -i(^n)aj:,{q, iuJn) + a_^{~q, -iu;„)a+g,(-^, -i'^n)) ■ 

Notice that Sp is an action for interacting bosons described by boson fields a(a+)^(g) satis- 
fying usual boson commutation relations [a^(^, '^^,(^)] = ^kk'^qq' [^a, ^/s] = [c^o, c^/j] = 
and with kinetic energies \k.q\/m and interaction term of form ~ v{q){a~^a+a'^a'^+aa+aa^). 
In this form the dynamics of the original fermion system is described completely in terms 
of boson fields (bosonized). Notice that we have so far restricted ourselves to the Gaussian 
approximation. Higher order interaction terms between bosons will appear in a cumulant 



expansion of the A^(g) fields We shall assume in the following that these high-order 
terms do not modify qualitatively the physics described by the Gaussian theory. Notice also 
that in the g — > limit, Ag(g) — 5(eg — and the usual "tomographic" bosonization 
procedure based on subdivision of Fermi surface into disjoint patches at small g*is recovered 
Our bosonization procedure can be viewed as a generalization of the tomographic 
bosonization method for small wave vector q to arbitrary values of g < A. 

To understand the nature of bosonization theory and Gaussian approximation, we first 
evaluate the free energy associated with Sp. Integrating out the Pj^iq) or a(a+)^(g) fields in 
Sp and using the fact that 

where [xo{q, p = icj^) and XoiQ^i^n) is the usual Lindhard function, we obtain 

Fp = ^ HXoM i^n)) + ^ E (1 - v{q)xo{q, iuJn)) , (H) 

k,q,iuJn g,i-^n 

where the first term in Fp is coming from the kinetic energy of the a(a+) bosons and the 
second term is the usual random-phase approximation (RPA) correction to free energy for 
interacting fermions. The presence of the RPA term in Fp suggests that our present Gaussian 
approximation is essentially the same as RPA theory for interacting fermions, except that 
the extra 'bosonization' assumption in our theory which gives rise to the extra kinetic term in 
Fp. The 'bosonized' action Sp assumes that the excitations in a fermion system can be fully 
represented by particle-hole pairs which are treated as independent bosons in the Gaussian 
apprioximation. Note that particle-hole pairs are not all independent in a fermion system 
because of Pauli exclusion principle. For example, two particle hole pairs f^fp and f^ff 
are not independent excitations because they both involve creation of fermions in state k. 
However, they are treated as independent bosons here as long a.sp ff. Notice, however that 
in one dimension the situation becomes different when the fermion spectrum is linearized 
near the Fermi surface. In this limit the entire particle-hole excitation spectrum can be 
represented rigorously by bosons p and the Gaussian approximation becomes "exact". At 
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higher dimensions, it is also beheved that the bosonization approximation is good as long 
as the interaction cutoff satisfies A « kp and as long as interaction with transverse 

gauge field is excluded 0. 

Despite these approximations in the Gaussian theory, the bosonized form of the action 
has the advantage that within the approximation the full excitation spectrum and the ground 
and excited states wavefunctions of the system can be obtained easily. This allows us to 
study the properties of the system in great detail, as we shall see in the following. 



III. PARTICLE-HOLE EXCITATION SPECTRUM 



The eigenstates and eigenvalue spectrum described by the action Sp can be obtained by 



diagonalizing the bosonized action (|ToD using a generalized Bogoliubov transformation. We 



introduce for each wave vector q the Bogoliubov transformation [|10 

k' 

k' 

and with correspondingly, 

k' 

7-ki-Q) = E ["k«-fc'(-^ - i^k'k^'U^ 



(12) 



(13) 



where we assume that the 'y{'y~^)j;{(f) operators diagonized the Hamiltonian, i.e. 

k 

where Ej:{q) are the eigen-energies and Eg is the ground-state energy of the system. Notice 
that a collective mode may appear in the system and is also included in the sum J^k- The 
matrix elements a and (3 satisfies the orthonormality condition 



E ['^kk"'^k'k" l^kk"(^k'k'' 
k" 

E ['^kk"f^k'k" ~ f^kk"'^k'k" 



^kit'^ 

0. 



(14) 



k" 



Writing down the equation of motions for ag(g)'s in terms of 7(7^)^(^'s |jTO|], we obtain the 
Bogohubov equations 

{EA^ - = "^11 ^(A^(^))^(A,., (g-))V|A,Kg-)A,.,(g-)|K.,„ + /5S„,^,), (15) 

k" 

it" 

Solving these equations we find that in general there exists two kinds of solutions: (i)particle- 
hole continuum, with Ej^{q) = |/c.g|/m, and (ii) collective modes, with energy Eoi^q) outside 
the particle-hole continuum satisfying the RPA eigenvalue equation l — v{q)xoi<l, Eo{q)) = 0. 
Most of the detailed mathematics can be found in Ref [ffOl. We obtain finally. 



^(A,^(g-))^(A,^,(g-))^|A,^(g-)A,„(g-)|^;e;;(g,|A:'.g1/m) 



kk' kk' T.d(\KA _ IM 

l^kk' 



-0(Ag(g-))^^(A,^,(g-))J|A,Kg-)A,^,(g-)|t;,^;(g,-|^'.g1/m) 



l^dlWjA _|_ \k-q\ \ 

for the particle-hole continuum spectrum /c', where Veff{q,uj) = — v{q)xo{(l,i^)) is 

the RPA screened interaction, and 



1 g(Ag(g1)v/|Ag(g-)| 
«a;o = . . . ..1 ' (16b) 



9a; 
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. 1 g(Ag(g-)).y|Ag(g-)| 

'?Xo('?.^) 



2 

uj=Eo{q) 



for the collective mode Eo{q). Notice that in the boson representation, the collective mode 
exists as an non-perturbative effect due to interaction and cannot be obtained from analytic 
continuation of the non-interacting boson modes. 

Next we examine the solutions of the bosonized Hamiltonian in the ^ oo limit. 

2 

First we consider the collective mode. Using the result that xo{q^^) ~^ the limit 

uj » kpg/m, where uq is the fermion density it is easy to see that in the limit oo, 
the collective mode frequencies are given by Eq^q) = cup, where up = (iziiioef.^i ^j^^ plasma 
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frequency. Notice that cup — oo as — oo, indicating that plasma oscillations are outside 
the physical spectrum in this limit. 

Despite the vanishing of collective excitation in the physical spectrum, the particle-hole 
excitation spectrum with excitation energies \k.q\/m survives in bosonization theory in the 
limit oo. In this limit Veff{q,\k'.q\lm) —l/xo{(l,W-'l\/^) the coefficients 

agg,'s and /?^g/'s remain regular, indicating that the particle- hole excitation spectrum is not 
qualitatively modified by effect of confinement. It is instructive to examine the charge fiuc- 
tuations carried by the particle-hole excitations by examining the commutator [p((f), 7^ (g*)]. 
In particular, the commutator will be zero if creating a particle-hole eigen-excitation does 
not introduce any charge fluctuations in the system, which is what we expect in the — 00 
limit. Using Eqs. @, (|^) and ( p.6a|) , and the usual boson commutation rules, it is straight- 
forward to show that 

[P(?1,7j(-g )] = S,,^0{A^M)^\A^,{q)\ X ] (17) 

1 - v{q)Xo{q, \k.q\/m) 

and vanishes in the limit ^ 00. 

Before ending this section let us examine our results obtained so far from bosonization 
theory. Within the Gaussian approximation, we flnd a RPA-like excitation spectrum with 
both collective modes and particle-hole excitations. As the coupling constant increases, 
the energy of the collective mode rises continuously to inflnity whereas the particle-hole 
excitation spectrum remains unaltered. The charge fluctuations carried by the particle-hole 
excitations are projected out gradually as increases, resulting in chargeless particle-hole 
excitations at g < A in the conflnement limit ^ 00. Notice that within the Gaussian 
approximation, the conflnement state analytically continues to the usual Fermi liquid state 
and there is no phase transition in between. The theory thus suggests that the conflnement 
state of a gas of fermions is a Fermi liquid state, with however c/iar^e/ess-quasi-particles 
constituting the Fermi liquid. It also suggests that this is a rather unusual Fermi liquid 
state, since bare fermions in the system carries charge e, and the quasi-particles must have 
vanishing overlap with bare fermions if they carry zero charge. 

9 



IV. SINGLE-PARTICLE PROPERTIES 



In usual bosonization theory for one-dimensional systems, the single particle properties of 
the system can be determined once a rigorous representation of the single-particle operator 
in terms of density operators Pl{<1) and pR{q) are obtained P]. In higher dimensions, this 
procedure becomes inadequate for two reasons: (l)The corresponding procedure requires 
that the bosonized representation of single-particle operator ipbi^) satisfies the commutation 
relations 



for all possible momenta k and q. We have not been able to find a representation which 
satisfies this criteria , and even if we can find such can representation, the theory is still 
approximate because in dimensions higher than one, the boson representation using Wigner 
operators is not exact and violates the Pauli exclusion principle, at least in Gaussian approx- 
imation. (2) more importantly, unlike in one dimension where the elementary excitations are 
collective density waves, we have seen in last section that in dimensions higher than one the 
particle-hole excitation spectrum is fermi-liquid like, implying that fermionic quasi-particles 
exist in dimensions higher than one. It is thus important to construct directly the quasi- 
particle operators in this case. 

To construct the quasi-particle operators we first consider the equation of motion of the 
bare fermion operator ip{f) = jd^k^~^^'^ fk at imaginary time, 

^ = ^^'^(^"') - ;^ E^(9)P(^)^^''>(^1 (18) 

q 



where 



{q)pA<i) = ll^\^Mn^km 

k 
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I ^ ^1 I ^ ^1 



(19) 



where v{q)pph{q) ~ 7^(^ + 7_a:(~9) describes the couphng of the particle-hole excitations 
to the fermion operator, and v{q)pc{q) ~ lo{(l) + 7o(~^) describes coupling of the collec- 
tive excitations (plasmons) to the fermion operator. In the — > oo limit, the interaction 
between fermions and particle-hole excitations are regular and finite, and the strong confine- 
ment effect shows up only in the interaction between fermions and collective mode excitations 
(plasmons). The effect of plasmons on the one-particle properties can be estimated perturba- 
tively by evaluating the plasmon contribution to fermion self-energy S and renormalization 
factor z to second order using Eqs. (|18D and (|19D. We find that S ~ Jl-i d'^qv{q), and 
Yjjijjp. Notice that both S and z goes to infinity as — > oo. Furthermore, the 
integrals carry also infrared divergence at dimensions d <2 for any finite e^. The divergence 
oi z aX d = 2 for finite was interpreted as signature of a marginal fermi liquid |12 . 



The divergence in single-particle self-energy indicates that the bare fermions is not a 
good starting point for constructing quasi-particle operators. To find a better starting point 
we first look for a canonical transformation for the single-particle operator which diagonalize 
the interaction term between fermions and plasmons. The kinetic energy of fermions and 
interaction with particle-hole excitations will be treated afterward. We obtain [|13| 



= e'^^^^V'(r), (20a) 

where 

E -T(7o^(?1-7o(-g1), (20b) 



ld/2 



q Eo{q) 



9X0 (g.^) 



du> 

which represents a fermion operator "dressed" by plasmon modes. It is straightforward to 
obtain the equation of motion of the "dressed" fermion ipq, 

= i (v'^Q(r) - - (V0(rl)2]^Q(f) - 2V0(f).V^Q(f)) (21) 

r 

where we have used the results [i?, 70(g)] = —EQ{q)jo{q) and [yj:{q),Jo{q)] = 0, etc. for 
k ^ and the usual commutator between bare fermions and density operators to derive 

11 



the above equation. We have also neglected constant energy terms coming from normal 
ordering of operators in Eq. (|2lD. It is clear from Eq. ( pT] ) that the direct coupling between 
fermions and plasmons is eliminated in the equation of motion of ipqlf). However, interaction 
between fermions and particle-hole excitations remains in the equation of motion. Moreover, 
an indirect coupling to plasmon is also generated from the " dressed" fermion kinetic energy 
term as is in the similar "small polaron" problem [l^. It is easy to see by direct power 



counting of in the 0(r) field that the coupling of the "dressed" fermion to plasmons 
through kinetic energy term is much weaker than the original fermion-plasmon coupling. 
In particular, we find that the self-energy correction of "dressed" fermions from (f){r) fields 
remains finite in the limit ^ oo. It is also straightforward to show that the infra-red 
divergence in fermion self-energy at two dimensions is removed for the dressed fermions 
because of the much weaker coupling to plasmons. 

These results suggest that the dressed fermion operators V'q(?^'s constitute a valid start- 
ing point to construct chargeless quasi-particles in the — > oo limit. The other interaction 
effects can be treated perturbatively. To show that this is indeed the case we first consider 
the commutation relation between charge and dressed fermion operators. It is straightfor- 
ward to show that 

m.i-Q(r')]=\ ^.yt'^if!^ (22) 



u)=Eo(q) 



du. 

which vanishes in the limit oo, when EQ{q) —>■ ujp —>■ oo, indicating that the dressed 

single particle operators '?/'Q(r)'s are indeed 'chargeless' in the — oo limit. To show that 
ipQ{r) and ipqir') defines independent quasi-particles when f ^ f' we check the commutation 
relation between the dressed fermion operators themselves. We obtain 

in the limits — oo and Ir* — -PI oo, where 



O 
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The vanishing of the commutator between different ipq operators separated by large distances 
at dimensions larger than one indicates that they can be used to construct independent quasi- 
particle operators when describing the dynamics of these systems at long distance. Notice 
that in one dimension, such a construction is not possible because of the long-rangeness of 
the commutation relation. In fact, the only fermionic operators which commute with density 
operators p(^'s are the ladder operators which raise or lower the number of particles in 
the system by one. There are only two independent ladder operators in the system, which 
change the number of left-going and right-going fermions and cannot be used to construct 
local quasi-particle excitations. 

Finally, it is also easy to show that 

< V'q(^V'q(^ >=< ip^{r)ip{r) > . 

The fact that the densities of bare fermions and dressed fermions are the same implys that 
the fermi surface volume of the dressed fermions is exactly the same as that of the bare 
fermions. Assuming that the dressed fermions form a free fermi gas, we also obtain 

< iiQ{r)ii{r) >~ e "0 -^^"^ " < il)^{r)il){r) >, 

to leading order in e^, which vanishes in the e — oo limit, indicating that the bare fermions 
and dressed fermions have zero wavefunction overlap, as is expected on physical ground. A 
similar calculation shows that the bare-fermion occupation number n{k) has no discontinuity 
across fermi surface in the — oo limit, in agreement with marginal fermi liquid picture. 

V. SUMMARY 

Using a bosonization approximation we studied in this paper a gas of fermions interacting 
via scalar potential v{q) = Ane'^/q'^ for q < A << kp- In particular we consider the e"^ ^ oo 
limit where the potential becomes confining. We note that because of the low momentum 
cutoff A << hp in our model, a crystal state is not expected to be formed and a new 



13 



treatment of the problem is necessary. Within a Gaussian approximation we find that the 
particle-hole excitation spectrum of the system is always fermi-liquid like, with the charge 
carried by the particle- hole excitation vanishing continuously in the e"^ oo limit. Based on 
this result chargeless fermionic operators are construct which, we believe can be used as the 
starting point for constructing quasi-particles in the system. The system can be considered 
as a 'marginal' fermi liquid where the chargeless quasi-particles have vanishing wavefunction 
overlap with bare fermions in the system. 

A major assumption we have made in our theory is that the Gaussian approximation 
for the Wigner bosons describes at least qualitatively correct the particle-hole excitation 
spectrum of the system. The Gaussian approximation is believed to be valid as long as 
the interaction potential is restriected to region of small momentum transfer q < A << kp 
0-^]. Notice that the approximation can be improved systematically via a cummulant 
expansion S. We find that to lowest order of the cummulant expansion there is no extra 



divergence introduced into the particle-hole excitation spectrum [jl4|. Notice also that we 
have considered here only the effect of longitudinal confining gauge fields in gas of fermions. 
Naively we expect that similar physics will be found in the presence of transverse gauge 
field. In the — oo limit, the chargeless quasi-particles we have constructed will decouple 
from the gauge fields and form a (marginal) fermi liquid. However, it is easy to show that 
the e*'^'-'"-* operator we have constructed is not able to "screen out" the transverse gauge field 
because it is constructed from density fiuctuations which couples to longitudinal gauge field 
only. The transverse gauge fields requires rather different treatment and we shall discuss it 
in a separate paper. 

The author thanks Prof. Zhao-bin Su for many helpful questions and comments. This 
work is supported by HKUGC through RGC grant HKUST6124/98P. 
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